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(54) METHODE POUR DETERMINER LA PERMEABIL1TE EQUIVALENTE D'UN RESEAU DE FRACTURE DANS UN 
w MILIEU SOUTERRAIN MULTI-COUCHES. 

(57) Methode de determination de Ea permeabilite de fis- 
sure equivalente d'un reseau de fissures dans un milieu 
souterrain multicouche a partir d'une representation 
connue de ce reseau. 

La permeabilite de fissure equivalente d'un reseau de fis- 
sures dans un milieu multicouche souterrain est determi- 
nee en discretisant. a I'aide d'une procedure specifique. 
chaque fissure (F) du reseau de fissures en elements de 
fissure (R) {tels que des rectangles par exemple) et en defi- 
nissant des noeuds N representant des elements de fis- 
sure interconnects dans chaque couche du milieu et de- 
terminant des ecoulements de fluides (par exemple des 
ecoulements en regime stable) a travers le reseau discre- 
tise en imposant des conditions de pression limites et des 
transmissivites de fluide a chaque couple de noeuds voi- 
sins. La methode permet de relier de maniere systemati- 
que des modeles de caracterisation de reservoir fissure a 
des simulateurs double porosite en vue de realiser une mo- 
detisation plus realiste d'une structure geologique souter- 
raine fissuree. 

Application par exemple en production petroliere par des 
ingenieurs de gisement en vue d'obtenir des predictions 
d'ecoulement fiables. 
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fluide mat nee- fissure regis par des forces capillaires. gravitationnelles. visqueuses ei des 
mecanismes compositionnels. amsi que des cchanges mairicc-mairicc ( simulateurs double 
permeabilite double porosiici. De numbreux exemples de techniques de lart anteneur 
som cites dans les references indiquee^ ci-apres. 

Thomas. L.K. et ai : "Fractured Reservoir Simulation". SPE Journal ( FeM ier 1 9£T, , 
42-54; 

Quandaile. P. et al : "Typical Features of a New Multipurpose Reservoir Simulator". 
SPE 16007. presence au 9eme Symposium sur la simulation de reservoir de San Antonio 
Texas. 1-4 revrier 19S7: 

Coats. K.H. ; "Implicit Compositional Simulation of Single-Porosity and Dual- 
Porosity Reservoirs". SPE 18427 presente au Symposium SPE sur la simulation de 
reservoir de Houston. Texas. 6-8 fevrier 1989. 

L'un des problemes que rencontrent les ingenieurs de gisemem consiste a doter ce 
modele de base de parametres en vue dobtenir des predictions d'ecoulement fiables. En 

15 particulier. les proprietes petrophysiques de base des fissures et des matrices, ainsi que la 
taille des blocs matriciels. doivent etre connues pour chaque cellule du simulateur 
d'ecoulement. Alors que la permeabilite matricielle peut etre evaluee a partir de caroties. 
la permeabilite du reseau de fissures contenu dans la cellule, cest-a-dire la permeabilite 
de fissure equivalence, ne peut pas etre evaluee simplemem et necessite la prise en 

20 compte de la geometrie et des proprietes du reseau de fissures. 

Une methode directe connue permet de determiner lecoulement en regime stable dans 
un reseau de fissures. Elle comprend M utilisation de mailles fines et regulieres 
conventionnelles discrepant a la fois les fissures et les blocs matriciels du volume de 
roche fissuree parallelepipedique considere. Pour differentes raisons, cette methode 
25 connue ne fournit pas de resultats fiables. sauf si le volume de roche fissuree est 
discretise en utilisant une maille comportant un nombre extremement eleve de cellules, ce 
qui necessite des moyens de calcul tres importants. 

Dautres modeles specifies servant a calculer les permeab.htes equ.valentes de 
reseaux fissures en 2D ou }0 sont egaiemcnt connus. par exempie par : 
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- Odl.ng. N.E. : "Pen^b.li.y of ■ Nalura| and SlmM 
and Tectonic Meddling and i[s AppliciKion to Pelmkum Cenk ^. \pp s ia| 
Publication 1 . 365-380. Elsevier. Norwegian- Petroleum Sucietv ,NPF, 



- Long. J.C.S. et al : "A Model for Steady Fluid Flow in Random Three-Dimensional 
Networks of Disc-Shaped Fractures '. Water Resources Research ,aou, vol ^ n 

S. I 105-1 115: 



- Cacas. M.C. et ul : 'Modeling Fracture Flow w, [h a Stochastic Discrete Fractun 
Network : Calibration and Validation. I. The Flow M ode r. Water Resources Researc. 
(mars 1990) vol.26, n"' 3: 



- B.llaux. D. : "Hydrogeologie des milieux fractures. Geometrie. connectivity et 
complement hydruulique". These de doctora. preencee a I'Eco.e Nationa.e Superieure 
des Mines de Paris, document du BRGM n° 186. Editions du BRGM. 1990: 

- Robinson. P.C. : "Connectivity. Row and Transport in networks Models of 
Fractured Media '. These de doctorac. St Catherine's College. Oxford University Ref ■ 
TP 1072. mai 1984. 

La presence invention concerne une methode de determination de la permeabi.ite de 
fissure equivalence dun reseau de fissures dans un.milieu multicouche souterratn. 

La methode est caracterisee en ce quelle comporte les etapes suivames : 

- discretisation du reseau de fissures en elemencs de fissure (te.s que des rectangles 
par exemple) et definition de noeuds represents les elements de fissure interconnects 

dans chaque couche du milieu, et 

- determination des ecoulements de n u ,des a travers le reseau d.scret.se en .mposan, 
des condu.ons de pression l.m.tes et des transmissivites de fluide a chaque coup.e de 

noeuds voisins. 

Plus prccisement. la methode est caracterisee en ce que: 
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- Ion divisc le milieu en une scrie de couches parallels se.cndam chacuno dans un 
plan de reference pcrpendiculaire a luxe de reference el dclimes chueune par une 
cooidonnee le longdudit axe. 

- on divise chaque fissure en une serie de rectangles limiies le long dudit axe de 
reference par deux couches adjaccntes ei on catalogue les rectangles en leur assoeiam 
des curacteres geometnques et physiques tels que coordonnees et tallies des reetancles et 
conductivites hydruuliques des fissures. 

- on positionne des noeuds dans chaque couche pour toutes les fissures 
interconnectees. et 

10 - pour tous les couples de noeuds vo.sins. on calcule des facteurs de transmissivite et 

resolution des equations decoulement en vue de determiner les permeubilites 
equivalentes du milieu dans trois directions orthogonales. 

Duns un mode de realisation prefere. la permeabilite equivalents du milieu comprend 
la determination directe du tenseur danisotropie de permeabilite equivalente et 
15 letalonnage de yaleurs absolues de permeabilite a partir de resultats d'essais de puits. 

La methode ainsi resumee permet de relier de maniere systematise des modeles de 
caracterisation de reservotrs fissures et des simuiateurs double porosite en vue de realiser 
une modelisation p!us realiste dune structure geologique souterraine fissuree. La 
methode peut etre mise en oeuvre par exemple en production petroliere par des 
20 ingenieurs de gisemem en vue dobtenir des predictions decoulement fiables. 

Dautres caructeristiques et avantages de la methode selon linvennon apparaTtront 
plus clairement a la lecture de la description ci-apres de modes de reahsanon donnes a 
titre dexemples non limitat.fs. en se referant aux schemas f.guram en annexe parmi 

lesquels : 

-J - la figure I monirc par exemple une image er. 3D dun rescau de fissures 

stochas, .quement obtenue a partir dobserva.ions e. de mesurcs prises sur un 

atflcurement de »jre>. 
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• leN figures :. 3 montrcn. une fissure divisee en une serie de rectangles R. 

- la figure 4 moi „re la ^ donn , cs ^ ^ ^ 

t issu rc. 

- la figure 5 momro !e mode prefer* de discretion dun p| a „ dc t , S sure. 

- les figures 6. 7. S « 9 illusiren, schematiquemen, le calcul des tae.eurs de 
,ransm,ssiv„e pour different positions de noeuds les uns par rapper, aux aucres ou avec 

une iimite. 

Le, permeability equivalen.es d un reseau de fissures 3D son, determine ci-apres e, 
u„l,sam une technique numerique basee sur la methode connue de "reseau de resistance" 
10 decn, par exemple dans ,e document d'Odling N.E. cite ci-avan, Dans ,a presente 
methode. on adme, que la motr.ee es, impermeable afin quel.e soi. compatible avec 
rapproche doubie porosi.e. Dans les simu.ateurs de reservo.r. les ecoulements matrice- 
fissure e, ma.rice-matrice son, en fait calcu.es separemen, par rapport aux elements 
au sein des fissures. 

"5 On admet que le reseau de fissures en 3D consider, represent, dans un volume e*a. a 

une ceUu.e de reservoi, h distri5ution ^ des ^ ^ ^ ^ 

caracteres de fissure du champ dans un mode.e de cassation. Le pnncipal objectif 
de ca.cu.3 decou.ements monophasiques dans le reseau de fissures en 3D est deva,uer 
.an.sotrop.e de permeabi.ite equivatente (Kv/Kh e, Ky/K.x) de ,a cellu.e de fissure 
cons.deree. qu, const.tue un parametre important regissant ,e componement des 
ecou,emen,s polyphasy dans les reservoirs. Les valeurs de permeabilue equ.va.ente 
obtenues par ces calcu.s pourraien, en pratique etre comparees aux resultats dessais de 
putts en vue d'eta.onner les caracteres de fissures tels que les conductivi.es nvdruuliques 
de fissure ,ou les ouvertures hvdrauliques equivalent, qu, ne peuven. etre de,ln,es que 
de maniere mediocre a priori. 

De plus, les resul.ats.de P ermeab,l,te equivalent* peuven. e.re unhses en vue V 
terminer un tenseur dc permeub„„e don, les pnncpa.c, directions permcen, u,, 
orientation opt.ma.e dc ,a maille du mode.e de reservo, Cependan, des cond.nons 
l-ites specfiques son, necessa.rcs pour ob.en.r ces informations. Des cond.nons linliles 



" 2757947 



sans eeoulement lateral imposes <ur les quatre faces lateral du volume 
parallciepipcd iq uc erud.e ne donncm pas acees aux termes non d.ugonaux du [cn»eur de 
nermeab.lue equivalent, aiors que do poientiels tou de> prcssion.s) a ^a^aln>n iineaire 
sui- Ics laces laterales permeuent d'imposer la direction du gradient de polemic! au scm 
5 du milieu ani.soirope et de deduire d.rectemem les termes de permeabiliie non dmaonaux 
Jos debits luteraux. 

Les techniques d'integrarion de donnees de fissuration naturelle dans des modelcs de 
reservoirs fissures sont bien connues. Les donnees de fissuration sont essentiellement de 
nature geometrique et comportant des mesures de densite. de longueur, dazimut et 
d'inciinaison des plans de fissure soit observes sur des affieurements. dans des calories de 
mines, sur des carottes, soit dedmts a pariir de diagraphies. Differents ensembles de 
fissures peuvent etre distingues et characterises par differentes distributions statistiques de 
leurs caracteres de fissure. Une fois les schemas de fissuration caracterises. on peut creer 
des reseaux numeriques de ces ensembles de fissures en utiiisant un precede stochastique 
respectant les distributions statistiques des parametres de fissures. De tels procedes sont 
decrits par cxemple dans les brevets FR-A-2.725.S 14: 2.725.794 ou 2 733 du 
demandeur. 
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La methode selon Invention sapplique a des images de structures geologiques 
flssurees presentant differentes tailles ou volumes et/ou situees en differents endroits. 
creees par un generateur de modele de fissure. La figure 1 presente une telle .maee. 

DONNEES D'ENTREE 

Avant de developper les procedures recommandees en vue de determiner les 
parametres hydrauliques equivalents d images de fissures en 3D. une etape important 
consiste a deTinir une structure commune de donnees demree pour ces images afm 
quelles puissent etre trauees independamment de fount de traitement utilise pour leur 
creation. 

Comme le monuent lex figures 2. .-. on admet que les fissures F son. xensiblemen, 
SC " icaleS '^-^n porpsndicubires aux l,m,.ex dc la coucho. Ceperulani. une ,nC-n,e 
S,rUCtUre di don ^ s &~ uppl.uuee aux Assures xocununc Icycrcmen. de la 
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xemcaie. Umaye 3D cm discrete- vemcalemcni en respectant la .irautlcai ion 
yeoloytque reelle lorsque eette information est dispon.bie. Dans le cas comra.re. on 
applique une discretisation arh.traire a ('image. Chaque couche honzontale L est 
caracterisee par sa coordonnee verticale zL dans le svs te me de reference de coordonnees 
«OX. OY. OZ). 



Une ser.e de rectangles R doit etre dermic pour chaque couche L. Chaque rectangle 
consiste en i.n element de plan de fissure compris entre les l.mites dune couche donnee. 
Ainsi. chaque fissure naturelle consist* en une scrie de rectangles superposes R ct on lui 
attribue une origine (origine de la fissure ). Chaque rectangle est defim par : 

i» - les trois coordonnees (xO. yO. zO> de longme du rectangle .0. Pour une fissure 

naturelle donnee, tous les po.nts d'origine des rectangles constitutes sont situes sur la 
me me ligne verticale (ou plus fort pendage) tracee depuis I'ongine de la fissure: 

- les coordonnees du vecteur d'unite horizontal i ( X H. yH) et du vecteur dumte 
verticai ] (xV\ y V) definissant I'orientation du rectangle dans le svsteme de reference de 

15 coordonnees. avec xVerticai et y Vertical egaux a zero en cas de ftssures vertices mais 
consideres comme des donnees d'entree afm de pouvoir traiter des fissures non 
verticales; 

- les deux longueurs horizontals algebriques I- et 1+ separant I'origine du rectangle et 
les deux limites laterales (venicales) de ce rectangle; 

:o - la hauteur h du rectangle, cest-a-dire la longueur du rectangle dans le sens ] . qui est 

lepaisseur de la couche si la discretisation en sens j correspond a la geolosie: 

- la conductive hydrauiique c dehvee de Implication de la loi de Darcy relative , 
lecoulemem dans les fissures (pour un gradient de pression ££. le debit au se,n de la 

ch .4 P 

fissure de hauteur h es, u etant la v.scosite du fluide,. La conductivity c est 

donnee P ar la relation c=k.u. ou k = a 2, , ; , cn ulilisanl la Kpri ^ mMi , m ^ 
fissures de Poiseu.llc, es, la permeabiiite mmnseque de la ii.w,re a son ,u,v,,„„c 
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h\draulique equivalente. La conductivity' hydraulique c e>t une \aleur de reference 
donnee pour une direction de la com mime horizontals maximale parallcle au >en.s de la 
fis>u re: 

- Ie> deux reciungieN -aipericur e! interieur voisin.s f 'R. LR: 

5 - la serie de fissures FS a lacueile appurticnt le rectangle considers; 

- Tangle d'orientation ct 0 de la direction de la contraime honzontale maximale prise a 
partir de I'axe (OX) dans le systeme de reference de coordonnees: 

- pour chaque serie de fissures, une table de correlation correlant I ) Tangle compris 
entre la direction de la contraime horizontale maximale et la direction de !a fissure 

!0 (azimut) avec 2) la conductivity hydraulique c ou i'ouverture hydraulique equivalente a 
precedemment definies. Les termes "horizontal" et "vertical" utilises dans ce contexte se 
rapportent aux directions respectivement paralleled et perpendiculaires aux limites de 
couches considerees ici comme etant horizontales. Les limices de couche discretisem les 
plans de fissure dans le sens "vertical". Fl convient de noter que les donnees d'entree 

15 menttonnees ci-dessus 1 ) conviennent pour tous les outils logiciels existants utilises pour 
caracteriser et generer des fissures et 2) qu'elles pourraient ete utiiisees pour discretiser 
un reseau de fissures legerement non verticales. cest-a-dire qui ne sont pas 
perpendiculaires aux limites de couche. 

MODES OPERA TO I RES 

20 Les modes operatoires et les tests de validation de la methode permettant de calculer 

I'anisotropie de permeabilite du reseau de fissures dans son ensemble sont deems ci-apres 
a partir de 1'image en 3D ainsi codifiee. La procedure numerique destinee a calculer les 
permeabilites equivalents d'un reseau de fissures 3D est egalement decrite. 

Le probleme consiste a trouver la distribution des debits dans le reseau pour les 
25 conditions from ie res suivuntes aux limites du volume paraMelepipedique etudie. c'e.st-a- 
dire le> pressions fixes imposees sur deux faces opposees et ies pressions variant de 
maniere Imeaire sur Ic.s quatre faces laterales (entre les valeurs imposees sur les deux 
autres faces). 
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Lcs principalis e tapes sont resume^s.ci-Jessous : 
I t Discretisation du rcseau 

En uhl.sam les definitions donnees pour la structure des donnees dentree. le reseat, de 
''• SSUre CSt ^^''^ NOUS b formc dune ^ "noeuds" N. chaque nocud clan, place 
UU Ctt " ,re de> ^'"onts dlntersecfon IS de deux rectangles R .ccs.-a-dire de deux plans 
de fissure au se,n dune couche donnee.. Comme le montre la figure 5. des noeuds 
supplemen.aires AN sont places au-dessus et au-dessous des noeuds precedents N , 
de represent dautres rectangles discretisant les fissures e. de minimis les longueur 
decoulement au .sein dune fissure donnee. Sur lu figure 5. BL est une Hmite lateral,- de 
deux cellules de fissure voisines. 

Une fois le reseau discrete, une procedure de tri est appliquee a ce reseat, de fissures 
en vue delimtner les noeuds isoles ou des groupes de noeuds sans lien avec l une des 
limi.es latera.es FL du volume 3D etud.e. car les fissures ainsi "tnees" ne comribuen, pas 
au transport du fluide et peuvent entraver les procedures de resolution mises en oeuvre 
en vue de trouver les pressions aux noeuds de fissures lors dun ecoulemen. en regime 
stable a travers le reseau. 

2) Caicul des transmissivites 

Un facteur de trunsmissivite T est calcule. pour chaque paire de noeuds associes en 
utilisant la relation : 

T _ ch _ Icah 
1 1 

ou c est la conductive hvdraulique de !a fissure, k la pcrmeabilW intnnseque de la 
fissure, a louverture de la fissure, h la hauteur de ia fissure et 1 la distance entre deux 
noeuds d'une fissure. 

DifTcren.es situation* do.vent etre envisages se.on la position respective des deux 
noeuds. Pour des noeuds con.enus au se,n de la meme couche ,,",,.6,. Ic taccur do 
transmissivite horizontal T est ob.enu d.rec.cment sous la forme dc la distance ,IUP, 
separan. les deux noeuds dans le sens de .ecoulement .fig.7, Pour JeN n „, Ut|s si(u , N 
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dans deux couches di.Veren.es .Fiy.^Je taceur do ,ra»>m,s>iw,e horizontal est la 
somme arithmetic des facteurs de transmissivite iT'+T" I relates aux deux elements du 
plan de fissure des cellules de fissure superposees. [I co,np„r,e unc k , n , Jlltfl „. 
decoulement egale a la demi-somme des deux epaisseurs de couche hi et h;. Pour les 
noeuds supplemcn.uires , c l, que defims prccedemment. relies par un seul plan de fissure, 
un seul tacteur de transmissivite est calcule pour cet element du plan de fissure. 

Le facteur de transmissivite T emre un noeud et une limite du volume 3D etudie est 
ex prime comme dans le cas de deu\ noeuds. selon les deux cas suivants. 

Pour une limite vemcale laterals. -le facteur de transmissivite T peut eire exprime 
directement pour un seul element du plan de Hssure (fig.8). e. sous la forme de la somme 
de deux transmissivites si deux plans de fissure relient le noeud et la limite. 

Pour une limite inferieure ou superieure horizontal le tacteur de transmissivite 
verticale petit etre exprime en considerant une longueur decoulement egale a la moitie de 
lepaisseur de la couche ( fig. 9). 

3) Equations d'ecoulemem 

En regime stable, un ecoulement monophasique incompressible a travers le reseau de 
fissure est determine en resoivant un serie de n equations, une pour chaque noeud. 
comme il est connu dans Part. Chaque equat.on exprime le fait que le debit total est nul a 
chaque noeud de fissure. Pour calculer un tenseur de permeabilite. on cons.dere qu'une 
pression constante est imposee a chacune des limites aval et amont. Une pression v ariant 
de maniere lineaire en fonct.on de la position entre les limites aval et amont est ,mposee. 

La matnce de permeabilite equivalente (Kijj determinee precedemment est 
diagonals pour calculer les directions principals de lecoulement avec les 
permeabilites equivalentes respect ives dans ces directions. 

En pratique, le problems se limite souvent a la decouverte des principales directions 
horizon, ales de lecoulement I' et V dans la mesure ou la direction P er P endtcula,re aux 
l.mi.es de couche , gencralemcn. verticals, est toujours prise comme I'axe des ,. Dans un 
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lei cos. seuls ,es .ennes ev,ra-dia W K, v et K vs doivent cue caleuK* ,ls peuvent 
e.re obtenus avec les conditions limitcs mixtes suivantes : 

- les ecoulements horizontal son, calculcs avec des .aces inKrieurcs et superior 
impermeahics e t des presses a variation lineaire sur les taces vertical parallel au 

> >en.s de lecoulemcm. 

- I ecoulemeni vertical est calcule avec routes les faces laterales impermeahles. 

Ainsi. un tenseur de permeab.Iue simple est obtenu. a purtir duqucl on P em 
facilemem dedu.re le> principals directions horizontals de 1 ecoulement U et V ; 

fK» Km 0 ^ 
0 0 K«J 
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Validation 



La methode a ete validee avec succes par rappon aux calculs decoulemem 
monophasique de reference ment.onnes c.-avam. realises avec un simulator de oisemen, 
conventionnel. Les calculs de reference ont ete obtenus sur de fines ma.Ues resul.eres 
discrepant les fissures, ainsi que les blocs matriciels du volume de roche fissuree 
paral.elepipedique consider*. Pour un sens decoulemem donne. une pression ^injection 
et une press.on de production fixes on. ete imposees sur les faces den.ree et de sortie et 
le debit resultant a ete calcule pour des conditions exemptes decoulemem lateral. 

On a realise trois etapes pour valider le calcul de : 

- la permeabilite venicale equivalent dun volume de roche traverse par une seule 
fissure, cette derniere etant representee par plusieurs noeuds correspond™, aux 
intersections avec de petites fissures isolees: 

- les permeab.lites honzontales equiv.len.es .dans une geometric d'ecoulement en 2D. 

et les principaux sens decoulemem; 
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- Ies permeab.Iites eqim alemes et ^isotropic de permeahiliie dans 
pre.seniani unu geometric J ecouleniem en .^D. 



un sen I reseau 



Les resultats obtenus pour la iro.sieme etape (pour une geomeme decoulemem en 
' D ' S ° m d ° nn ^ Jans ,tf t:ib,eaLt "-ap«N. Pour Ies directions horisontalcs decoulemem. 
une .solution anaiynque de reference pcut egalemem etre caleulee puisque la geometric 
de lecoulement esr un ecoulemeni :D dans ces directions < !u L'eometne decoulemem en 

concerne la di reel ion z>. 



Permeabilites 
equivalentes ( md) 


Simulation 
MAILLE FINE 


METHODE 
PRESENTEE 


Solution 
ANAL YTIQUE 


Kx 


0.1 19 


0.120 


0.120 


Ky 


0.224 


0.227 


0.226 ' 


Kz 

I 


0 2^ 
— 


0.267 




Anisotropic 
Kz/(KxKy)0-5 


1,56 


1.62 





II apparait clairement que Ies resultats obtenus avec la methode presentee sonr tres 
proches des valeurs correspondantes obienues avec la soiut.on analytique et | a simulation 
a maille fine pour Ies directions X et Y. 

De plus, la difference entre Ies valeurs de permeabilite equivalente venicale 
concernant lecoulement 3D reste acceptable. Ainsi. le rapport danisotropie. egat a 1.6. 
est predit de maniere satisfa.sante par la methode avec un nombre tres limite de cellules. 

La methode selon Invention, qui fourn.t une representation facilemem transposable 
dun reseau de fissures nuturel. se prete bien au calcul decoulemcnts au sen de fissures. 
Elle peut egalement etre utile pour ameliorer Hmage dorigine du reseau de fissures. Une 
telle image est en fait obtenue a panir dun generateur de fissure stoehastique uti | isam 
comme donnees d entree Ies resultats de Integration de donnees de fissuration urchnees 
dans un modele de caraclcnsahon de llssure comme ,1 est decnt dans Ies brevets FR-A- 
17;5.SI4: ;.725.794 , 1U 2. 7^3.07? du den>andeur deja mentionnos. De .elles iniaucv 
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une tois discret.sees par la methode selon (invention, peuve.u Cm re facile meni mod. flees 
aim d'etre adaptees aux regies »eolojiiques. Par e.xemple. (interrupt ion >wemar N ue 
d une fissure donnee centre une autre sene de fissures peut etre prise en eompte duns 
I'ima«e d'ongine en annufant les elements du plan de fissure dune sene donnee > etendant 
au-dela des fissures entrecoupees de lautre sene. 
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REVINDICATIONS 

I ) Methode permeitant de determiner la permeabilite do fissure equivaiente d un 
reseau de fissures dans un milieu multicouche sou terrain fissure, a p- M in d une 
representation predeierminee dudit reseau. comportant les etape^ de : 

- discretisation de chaque fissure <F) du reseau de fissures en elements de fissure et de 
definition de noeuds (Nt represemaru des elements de fissure interconnected dans chaque 
couche du milieu, et 

- determination d'ecouiemems de tluides a travers le reseau de fissures en imposam 
des conditions limites de pression. et de transmissivites de fluides pour chaque couple de 
noeuds voisins. 

2) Methode selon la revendication I, caracterisee en ce que: 

- I'on divise !e milieu en une serie de couches parallels setendant chacune dans un 
plan de reference perpendiculaire a i'axe de reference et defmies chacune par une 
coordonnee le long dudit axe. 

- on divise chaque fissure en une serie de rectangles limites le long dudit axe de 
reference par deux couches adjacentes et on catalogue les rectangles en leur associant 
des caracteres geometriques et physiques tels que coordonnees et tailles des rectangles et 
conductiv ites hydrauliques des fissures, 

- on positionne- des noeuds dans chaque couche pour toutes les fissures 
interconnected, et 

- pour tous les couples de noeuds voisins. on calcule des facteurs de transmissivite et 
resolution des equations decoulement en vue de determiner les permeabilites 
equivalentes du milieu dans trois directions orthogonales. 

3) Methode selon i'une des revindications precedemes. caracterisee en ce que la 
permeabilite equivaiente du milieu comprend ta determination directe dun teiiseur de 
permcabiliie equivalent et letalonna-e de valeurs absolues de permeabilite a partir dc 
resultuts d'es.sais de puits. 
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4) Meihode selon Tune des revendicauons precedences, earacterisce en ce que lesdiis 
caracteres geomeiriques et physiques soni les coordonnees et les dimensions des 
rectangles ei les conducm iies hydrauliques des fissures. 
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FIG.4 
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